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The QCD hardness scale for secondary particles (h) production in semi- 
inclusive deep inelastic scattering (DIS), ep — ► e'Xh, gradually decreases from 
Q 2 , the photon virtuality which determines the hard scale in the virtual pho- 

\Q ■ ton (current) fragmentation region to a soft, hadronic, scale in the proton 

fragmentation region. This suggests similarity of the inclusive spectra of 

q^ . leading protons and neutrons, h = p,n, in high energy hadron-proton and 



virtual photon-proton collisions. We explore this similarity extending to the 

Mh. DIS regime the nonperturbative peripheral mechanisms of inelastic scattering 

traditionally used in hadronic interactions to explain fast nucleons production. 



^ . While the production of leading neutrons is known to be exhausted by DIS 

off charged pions, the production of leading protons by DIS off neutral pions 
must be supplemented by a substantial contribution from isoscalar reggeon 
(/o) exchange extrapolated down to moderate values of xl- We comment 
on the x and Q 2 dependence of leading proton production as a probe of a 
universal pattern of the x, Q 2 evolution of the nucleon and meson (reggeon) 
structure functions at small x. 



In deep inelastic ep scattering, according to the standard QCD description of hadroniza- 
tion, the proton fragmentation region is very different from the current and/or photon frag- 
mentation region. Namely, the virtuality of partons in generalized ladder diagrams gradually 
decreases from the hard scale Q 2 of the struck parton in the current fragmentation region 
to the soft, hadronic, scale of the parent parton in the proton fragmentation region. 

Although until quite recently experimental data on proton fragmentation were scarce 
Q, presently the ZEUS and HI leading proton spectrometers (LPS) and forward neutron 
calorimeters (FNC) are operational and are amassing data on leading proton and neutron 
production |||J. 

Whereas popular Monte-Carlo implementations of perturbative QCD (Ariadne ||, Her- 
wig || and others) are very successful in the photon fragmentation region (for a recent 
review see 0) , a purely perturbative description of the proton fragmentation region is not 
yet possible and the current versions of Monte Carlo hadronization models underestimate 
the yield of fast secondary nucleons || . Traditionally, leading proton production in inelastic 
collisions is modeled via nonperturbative peripheral interactions. Such peripheral models 
were quite successful when applied to hadronic collisions (for a review see ) . A well known 
example of this type of processes is the interaction of projectiles with pions from the chiral 
mesonic substructure of the proton. We recall that in hadronic reactions, the pion exchange 
mechanism with absorption exhausts the cross section for the leading neutron production 



in pp — > Xn, 7rn — > Xp, pn — > Xn (see |H| and references therein). There is at present a 



mounting evidence of the importance of the chiral mesonic substructure of the nucleon also 



in DIS (see |TT]] and the reference therein). Here let us mention a successful explanation of 



1 Before the ZEUS/HI LPS era, leading protons in DIS have been studied only in the fixed-target 
bubble-chamber neutrino experiments [Q]. Only large values of the Bjorken x were accessible, 
leading to a strong kinematical bias in the leading proton spectra ||. The small- x data from 
HERA are free of this kinematical bias. 



the Gottfried Sum Rule violation (see for instance JT2]| and references therein ) and of the 



d-u asymmetry as seen in pp and pn Drell-Yan production [|nj. Exactly the same physics 
opens an exciting possibility to study the pion structure function at HERA. Here one can 
separate DIS off nearly on-mass-shell tt + by triggering on fast neutrons from semi-inclusive 
ep — > e'Xn reaction ||TB[ (see also ]TT|-|TTS|] ) . In this communication we explore to which 
extent leading proton production in DIS can be understood quantitatively within peripheral 
mechanisms. 

Following the conventions for diffractive DIS [pOJI , we define the semi-inclusive structure 
function 

; %^^ [2 - 2 ^r^ , '^ i - M2) - (1) 

Here a is the electro-magnetic fine structure constant, z is the fraction of the light-cone 
momentum of the beam proton carried by the outgoing proton z = ^+- (the same quan- 
tity is often denoted by xl), t is the (p,p') four- momentum transfer squared, x, y, Q 2 and 
(3 = M 2 +0 2 — i£j are the standard DIS variables and R = a L /a T . In the present analysis 
we focus on leading protons with 0.6 < z < 0.9. In our extension of the peripheral models 
of fragmentation used in hadronic reactions to lepton DIS at HERA, we consider four mech- 
anisms of leading proton production (Fig. 1): a) diffractive production of protons (pomeron 
IP exchange) which dominates at z — > 1, and constitutes a background to fragmentation in 
non-diffractive DIS at z < 0.9 — 0.95; b) spectator protons from the fragmentation of the nN 
Fock state of the physical proton produced by DIS off virtual 7r° (pion-exchange mechanism); 
c) protons produced as decay products of fast baryon resonances of which the A production 
via pion-exchange is a typical, and predominant, source; d) reggeized heavy meson (reggeon 
R) exchange contribution (predominantly the isoscalar reggeon, R = f , exchange). A pre- 
liminary evaluation of the first two mechanisms has been done in [jnj. The background 



2 For related work on electro-magnetic properties of nucleons and TF-boson and jet production in 



nucleon-nucleon collisions see [13,14|. 



from pion and reggeon exchange to the dominant pomeron exchange at z > 0.95 has been 
discussed recently in [p0|-|23|j, and we partially use the results of these works. Diffraction 



excitation of the proton into high-mass states also contributes to leading proton production, 
and we shall comment on this small contribution following the considerations in ||24| . 

Under approximations to be specified below, the contributions of all four mecha- 
nisms to the semi-inclusive structure function can be written in the factorized form 
(i = IP,7r p,7rA,/ ): 

F 2 (4) 0M,/3,Q 2 ) =$> 2 (4) (^,t,/3,Q 2 ) =J2fi(z,t) ■ F*((3,Q 2 ) , (2) 

i i 

where F%(/3, Q 2 ) is the structure function of the exchanged object (pion, pomeron, reggeon), 
fi(z,t) is its flux factor and (3 is the Bjorken variable for DIS off the exchanged object. 

We start our discussion with the pion exchange mechanism. In this case F^fl, Q 2 ) is the 
structure function of the physical pion and the flux factor is given by 

W""'ff"-" Ty ' <3) 

Strictly speaking, Eq. (|3]) holds in the plane wave impulse approximation. A recent analysis 

|T0| has shown that absorption corrections to pion exchange in DIS are small and can be 

neglected for the purposes of the present analysis. Also, the off-mass shell extrapolation 

effects are marginal and the on-mass shell pion structure function can be used. Important 

consistency check is provided by the simultaneous description of the hadronic leading nucleon 

data. The results for DIS in the interesting region of 0.6 < z < 0.9 only marginally depend 

on whether the light-cone or Regge parameterization of \F n ^(z, t)\ 2 are used (for a detailed 

discussion concerning the choice of the form factor see Refs. fl0| , |l2"|| ). 

Production of fast A's is also known to be dominated by pion exchange. For A ++ 

production the flux factor is given by 

, / .x _ ^9 P n-A++ (1 „ A(m A + m N ) 2 - t) 2 ((m A - m N ) 2 - t)\F nA (z,t)\ 2 
16tt 2 6m 2 N m A (t — m^) 2 

Contributions from A + and A production can be included using the familiar isospin re- 
lations (see for instance |12 , 24[ ). In the simplest one-pion exchange approximation, the 



polarization state of the produced A's is such that the A — > ttN the decay angular distri- 
bution in the Gottfried- Jackson (t-channel) frame equals: 



w(9j, <Pty) = 1/4 • (1 + 3cos 2 6j) ■ Y m (<p TY ) 



where 8j and and 4>ty are the so-called Jackson and Treiman-Yang angles, respectively. 
Absorptive correction modify slightly this simple form | 25fl , but these corrections can be 



neglected for the purposes of the present analysis, since both the z— and t— spectra of 
decay protons only weakly depend on the decay angular distributions. 

For the diffractive e + p — > e' + p' + X reaction, our semi-inclusive structure 
function coincides with the pomeron component of the diffractive structure function, 
F^\jPN;z,t,/3,Q 2 ) = F 2 D ^\z,t,/3,Q 2 ). At z < 0.9, diffractive DIS is a small background 
to non- diffractive DIS and a somewhat simplified description is justified. Since the ZEUS 
data have |3j x < 10~ 3 and z < 0.9 then /3 is quite small, (3 < 2 • 10 -3 — 10 -2 and it has 
been argued |2T)||2"T| ] that at such a small j3 one expects the factorization 



F 2 Di4 \z,t,(3,Q 2 ) = Mz,t) ■ F*(f3,Q 2 ) . (5) 

The normalization of the pomeron flux factor fjp(z,t) and the pomeron structure function 
F 2 IP (/3, Q 2 ) is a matter of convention, and only the product of the two is well defined. To be 
specific, we use the triple-Regge parameterization for the flux factor 

M*,t) = 87r 2 (1 1 _ z) ( 1 -^) 2(1 " aip(t))G! ip( t )» ( 6 ) 



where Gjp(t) = G exp(Bjpt) with G = 21.2 mb J2"T|-[2"3| from the Regge decomposition 



of the NN total cross sections |28j and Bjp = 3.8GeV -2 according to the triple-Regge 
analysis of hadronic diffraction scattering [^,^9p|3l[|. For z < 0.9, the specific Regge effects 
coming from (1 — z) 2 ^ aip ^ and from the ^-dependence of the pomeron trajectory ajp(t) 
are marginal and, besides the standard factor j^-, the main z dependence of the flux comes 

from the kinematical boundary \t\ > \t\ m i n = — in the form factor G(t). In principle 

F 2 -^> (z,t, P,Q 2 ) can be derived from experimental data on diffractive DIS, but currently 



for (3 % 2 • 10 — 10 the pomeron structure function stays basically unknown. It has 



been argued, |26| , that at small (3 the conventional DGLAP evolution holds for the pomeron 
structure function giving a /3, and Q 2 dependence of Ff(/3, Q 2 ) similar to that of F£((3, Q 2 ) 
(see for instance |S7|]). On the other hand, the triple-pomeron formula with soft pomerons 
gives the scaling prediction F^(j3,Q 2 ) = Cjp(3~ om (the normalization Cip = 0.026 has 
been fitted |23|| to the HI experimental data pOfl). We use these two models to check the 



sensitivity of the leading proton spectra to the evolution in (3 and Q 2 . 

The reggeon exchange is an important ingredient of the triple- Regge phenomenology of 
hadronic diffraction, although its strength is not very well known P, 29], 30] . The triple- Regge 
parameterization for the reggeon flux is 

fn(z,t) = ^- 2 (l-z) 1 - 2a ^G R (t) , (7) 

where or(£) is the reggeon trajectory. We take Gr(£) = Gn(0)-exp(B R t), where for the dom- 
inant /o _exc hange Gr(0) = 76 mb |H],^] and Br = 4 GeV -2 , which is consistent with the 
data on leading proton production in pp collisions |£J . Triple- Regge considerations in con- 
junction with fits to the NN total cross sections [^J suggest the isovector reggeon exchange 
to be much weaker than the isoscalar /o exchange |2"T| , |2"B|] . The reggeon structure function 
is basically unknown. The extension of microscopic analysis |26| to reggeons suggests that 



the j3 and Q 2 dependence of this structure function at small (3 must be similar to that of 
the pion structure function, F^((3, Q 2 ) ~ F£((3, Q 2 ) and/or the pomeron structure function, 
F^(j3,Q 2 ) ~ Ff{f3,Q 2 ). Arguably the gross features of the small-/?, large-Q 2 behavior of 
F^ifl, Q 2 ) should be similar for all targets i. If the extrapolation along the Regge trajectory 
to the particle pole t = m 2 r were possible, one could have related F^(f3, Q 2 ) to the /o meson 
structure function, which at small f3 is expected to be similar to i^(/3, Q 2 )- Going from the 
particle pole t = ml to the scattering region t < brings the off-mass shell suppression 
in, and it is natural to expect -F 2 H (/3, Q 2 ) < F|"(/3,Q 2 ). The triple-Regge phenomenology of 
hadronic diffraction suggests the suppression factor A/ ~ 0.5 with a large uncertainty pi 



(because of different normalization of the flux in |ZT| and [ 23j , the estimate of A/ in |21| must 



be taken with the factor |). Similarly to the pomeron structure function, the triple- Regge 
formalism with soft pomerons gives the scaling F^(j3, Q 2 ) = Cr/3~ ' 08 . Eventually, with high 
precision data on diffractive DIS, one would be able to evaluate F^(j3, Q 2 ) directly from the 
reggeon background to pomeron exchange at z > 0.95. 

The single particle inclusive (z, t)-spectrum of protons is defined as R(z,t,x,Q 2 ) = 
F^'(z,t, P,Q 2 )/F 2p (x,Q 2 ). A fully differential study of R(z,t,x,Q 2 ) is not yet possible 
with the limited statistics of the preliminary ZEUS data |§. The data were collected within 
the following experimental cuts fl exp : 0.6 < z < 0.9, \t\ min < \t\ < 0.5GeV 2 , 10~ 4 < x < 10~ 3 
and 4 < Q 2 < Q 2 max , where Q 2 max is the maximal kinematically attainable Q 2 . Within these 
cuts the fraction of events with leading proton is given by: 

p _ V^ pi _ y^ ylexp) ( o\ 

exp *? exp V ^"(n-p) ' [) 

where 

fZmax ftmax fXmax (Qmax n OXT 

Aa\Q exp ) = dz dt dx I dQ 2 , (9) 

J z min Jtmin ** x min ^ min tlJUtlVcg iXZiXb 

f x max rQ 2 max n d(J t0t 

Aa tot (n exp ) = dx I dQ 2 —— , 10 

Jx mm JQ 2 min dxdQ 2 

and the subscript i stands for one of the mechanisms shown in Fig.l. 

As emphasized above, the pion, pomeron and reggeon structure functions are unknown 
in the j3 region considered in our present analysis. For a reference evaluation of R l exp , we 
take the GRV parameterization for the ((3, Q 2 ) evolution of the pion structure function 
the flux of pions evaluated in the light-cone model for the chiral structure of the 



nucleon [|l!| (the Regge parameterization leads to very similar result ||10|| ) and the triple- 
Regge model parameterization F^(l3,Q 2 ) = 0.026/3 -008 described above. For the proton 
structure function, which enters the evaluation of the denominator (|H]) of the ratio (|8|), one 
can use any convenient fit to the HERA data. In the present analysis we take the GRV 
parameterization [Q. The practical calculations have been performed with a Monte Carlo 
implementation of the above formalism. As a result of our analysis, we find R n o p (Q exp ) = 2%, 



RwA^exp) = 0.9% and the tail of the pomeron exchange contribution p2| gives Rjp p (Q e xp) 
= 1.2%, so that -Ri+2+3(^ex P ) = 4.05%. Note that QCD hadronization models (Ariadne, 
Herwig) were never meant to describe the nonperturbative proton fragmentation region; for 
example Ariadne || gives R(Q exp ) in the per mill range and a similar under-prediction for 
the production of leading neutrons ||. ^From the comparison with the ZEUS experimental 
result, Rex P US = 9.2 ± 1.7 %(stat. only) ||, we conclude that about 5 % of the missing 
strength must be attributed to the reggeon exchange. In the triple-Regge scaling model, 
F^(/3, Q 2 ) = Cr/3~ ' 08 , this requires Cr = 0.12 within a factor of 1.5 uncertainty. 

The importance of different mechanisms can be better seen from the z-dependence of the 
ratio R exp (z) defined for the experimental (t,x,Q 2 ) range as shown in Fig. 2. Clearly the 
importance of the reggeon exchange can be seen from the figure. With the set of parameters 
specified above, the reggeon contribution makes R exp (z) approximately flat at z < 0.9, in 



close similarity to a flat ^-spectrum of leading protons in hadronic interactions ||32|| . The 
preliminary HI results are also consistent with the flat z-spectrum ||. 

The z-spectrum of leading nucleons from diffraction double dissociation (DD) has been 
studied in ]23||. It can be isolated experimentally by the rapidity gap (GAPCUT) selection 



method ||. An extension of the analysis J23J] to leading protons shows that ~ 70% of DD 
events have leading protons, mainly produced by excitations of the Nttit and high mass 
continuum states. Roughly ~ 50% of final state protons have z > 0.6. Since DD constitutes 
~ 2% of the DIS events, only ~ 0.7% have a leading proton with z > 0.6 generated by this 
mechanism. DD is therefore a small, /(GAPCUT) = RoDiS^exp) / ' Rexp ~ 7%, background 
to the dominant non- diffract ive production mechanism, in good agreement with the ZEUS 



findings. The LEPT06.5 'soft color interaction' model |35| which, unlike Ariadne and/or 
Herwig, is supposed to describe all aspects of DIS including leading proton production, 
over-predicts the fraction of GAPCUT events: /(GAPCUT)^20-30%. The observed leading 
proton ^-distribution of the GAPCUT sample [[J is also consistent with the z spectrum of 
protons generated in proton dissociation into Nitit and continuum states as shown in Fig. 2 



of [pijl - R can readily be included in the analysis of higher precision data. 

8 



This evaluation of the reggeon exchange from fragmentation into protons at z < 0.9 is 
consistent within a factor of 2 with estimates of the reggeon background to pomeron exchange 
in the diffractive region of z > 0.95 [p^| . A caveat in comparing these two extreme regions is 
the possible reggeon-pomeron interference contribution ~ i— , which can be substantial in 
the diffractive domain and small in the fragmentation region z < 0.9. A combined analysis of 
high precision fragmentation and diffractive data would be the best way to fix the reggeon- 
pomeron interference contribution, but such an involved phenomenology is not warranted 
with the presently available data. Note also that the ZEUS data are preliminary and lacking 
the evaluation of systematic errors. 

In Fig. 3 we show the slope b(z) of the ^-distributions defined in the experimental range of 
(x, Q 2 ) (R(z, t) oc exp(b(z)t)) . The slope of the reggeon trajectory is large, a' R = 0.9 GeV -2 , 
and for pure reggeon exchange contribution quite a substantial rise of the slope is expected: 
bii(z) = Br + 2a' R log tz-. Similar growth of the slope is expected also for the pion exchange 
contribution. The increase of the slope at large z is tamed by the small diffraction slope 
of the pomeron contribution. The parameter B R is poorly known and the leading proton 
spectrum offers the best possibility for its determination. Our results for b(z) obtained with 
Br = 4 GeV~ 2 are close to the slope of the ^-dependence for leading protons observed in pp 
collisions ||32|| . 

The (x, Q 2 )-dependence of different mechanisms is controlled by the ratios pi(x,Q 2 ) = 
F-2(i~;iQ 2 )/F2p(. x iQ 2 )- m the scenario with the scaling soft pomeron/reggeon structure 
functions, Pjp,r(x, Q 2 ) decreases with rising Q 2 and/or decreasing x, because of the scaling 
violations and steep ^-dependence of the proton structure function F 2p (x, Q 2 ). Fig. 4 shows 
that in such a scaling scenario (no QCD evolution for F 2 ' (/3,Q 2 )), one would expect 
significant dependence of the leading proton production on both x and Q 2 . On the other 
hand, if Fjp jR (P,Q 2 ) satisfies the conventional (P,Q 2 ) evolution at small (3 PS] , one would 



expect very weak (x, Q 2 ) dependence of the leading proton spectrum. This stays true also 
in the real photoproduction limit. In the present analysis we model the evolution effects 
by taking F 2 ' (/3, Q 2 ) = Aip i r/|"(/3, Q 2 ) with the GRV pion structure function. We adjust 

9 



Aip = 0.2 and A/ = 0.5 as it was evaluated in |^TJ so that we reproduce the same R l as 
with the scaling (no evolution) scenario within the ZEUS kinematical cuts. 

In the conventional evolution scenario we indeed find a very weak x and Q 2 dependence of 
the leading proton spectra. The preliminary ZEUS data || better agree with this scenario, 
although the error bars are still rather large. The preliminary HI data || on the ^-integrated 
F 2 (z, t, x, Q 2 ) also support the conventional evolution scenario. 

We conclude that the salient features of fragmentation into leading protons can be un- 
derstood quantitatively in terms of peripheral mechanisms extended to the DIS regime. The 
experimentally observed similarity of the leading proton spectra in pp collisions and ep DIS 
is a natural consequence of these mechanisms. We emphasize that our approach has the 
capability of a unified description of diffractive DIS at z > 0.9 and of fragmentation into 
protons in non-diffractive DIS. Of the four sources of leading protons pion exchange can 
be experimentally determined using neutron tagged DIS. Experimental confirmation of our 
estimate for this process will lend strong support also for our evaluation of the A contri- 
bution. The pomeron exchange background can be inferred from diffractive DIS. Finally, 
the reggeon exchange mechanism of fragmentation can also be tested in diffractive DIS. The 
combined analysis of the high precision leading proton data and diffractive DIS data makes 
possible a determination of the reggeon-pomeron interference effects, which has not been 
accomplished with the hadronic diffraction data P, p9~|j30"l| . 

A comparison of the soft pomeron no-evolution (unrealistic though it is) and conventional 
evolution scenarios for the reggeon structure function shows that the high precision leading 
proton spectrum offers an interesting test of the universality of the QCD evolution properties 
of structure functions at small x (For a related discussion of the pion exchange mechanism 



within Veneziano's fracture function [[I(J context see P7|.). 

We conclude with the comment that similar fragmentation mechanisms may be at 
work also at smaller z, where DIS on the multi-pion Fock states of the physical nucleon, 
(n,7r)N, (n7r)A, (wr)iV"*, (wr)A*, provides a natural mechanism for slowing down secondary 
protons. In the spirit of the above discussion, the weak x, Q 2 dependence must hold also 

10 



for slower protons. Similar arguments hold for the fragmentation of protons into hyperons 

(A, £,...). 
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Figure captions 

Fig.l 
Peripheral mechanisms of leading proton production. 

Fig.2 

The fraction (in per cent) of DIS events with a leading proton in a given z bin (Az = 0.03) 
predicted by our model (thick solid curve) in comparison with the ZEUS preliminary data 
H]. The contributions of four mechanisms of Fig. 1 are shown separately: the thin solid line 
shows the pomeron-exchange contribution, the long-dashed curve is for the pion-exchange 
contribution, the dashed curve shows protons from the A production and the dotted curve 
is for the reggeon-exchange component. 

Fig.3 

The slope of the t-distributions predicted by the model is compared with the ZEUS prelim- 
inary data 0. 

Fig.4a 

The sensitivity of the fraction of DIS events containing leading protons within ZEUS cuts to 
the Q 2 evolution effects for the two different scenarios for the pomeron and reggeon structure 
function: the curves marked by filled circles are for the no-evolution soft pomeron model, 
the unmarked curves show the results for the conventional QCD evolution scenario modeled 
by the GRV pion structure function. The legend of curves is the same as in Fig. 1. 

Fig.4b 

The same as Fig. 4a, but for the ^-dependence of the fraction of DIS events containing 

leading protons within ZEUS cuts. 
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